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ABSTRACT 


GRAVITATIONAL HARMONICS FROM SHALLOW RESONANT ORBITS 


C.A. Wagner 
S.M. Klosko 


Until very recently, there has been no identification 
of the significant gravitational constraints on the many 
common orbits in shallow resonance. Without them it is 
difficult to compare results derived for different sets of 
harmonics from different orbits. With them it is possible 
to extend these results to any degree without reintegration 
of the orbits. Five such (strong) constraints have been derived 
for the GEOS II orbit (order 13, to 30th degree) whose princi- 
pal resonant period is 6 days. The constraints explain the 
sinusoidal variation with argument of perigee of a lumped 
harmonic found from 41 6-day arcs of optical and laser data 
in 1968-69. For example, the constant terms derived are: 


10^(38.1,-55.9) = -.872(C,S)^3 ^3 


(C,S)i5^13 + 


in terms of fully normalized spherical harmonics. 


The condition equations, derived from elementary pertur- 
bation theory are shown to account for almost all (>98%) of 
the resonant information in the tracking data. They agree 
well with recent gravitational models wliich include sub- 
stantial amounts of GEOS II tracking data. 
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INTRODUCTION 


We can calculate the amplituder and phases of the geo- 
potential perturbations from the formulas given by Kaula, 1966. 
He has expressed the potential entirely in terms of Kepler 
elements as 


Y = t. Z t Z ^"o«,rwi Where 
It m p q 


a 


^impq ~ — 


a 


e 

i+r 




cos 

sin 


C£-m) even 
C£-m) odd Cl) 


and 


*'£mpq "" C^>-2p)w + (£-2p+q]M+raCSi-0-XJlm) . 


In Equation (1) ; 

a,e,l,a) s'i, and M are the mean Kepler elements; 

semi-major axis, eccentricity, inclination, 

« 

argument of perigee, right ascension of the 
ascending node and mean anomaly respectively; 

a^ is the semi-major axis of the earth; 

u is the gravitational constant of the earth; 

0 is the rotation rate of the earth; and 

mXp is tan"^ — ^ . 

•"Jim 


"T 



The indexes it,m are the degree and order of a fully nrnnalized 
spherical harmonic term; the p,q quantities identify a par- 
ticular component of tliat term. The function (fully 

normalized) and arise when the potential is converted 

from position coordinates to Kepler elements. The 
functions are identical to Hansen's coefficients. The 
functions are sinusoidal with wave length about 2tt/ (Jl-m+1) . 

In terms of Kepler elements, the equations of satellite 
motion are [Kaula, 1966, p, 29] 

dt na 3M 

^ H _ (i> B«) t/z ^ 

dt na*e 3M na*e 3« 

^ ^ ^ cos i 3V ^ (1-e) 3V 

dt na* (1-e*) i 31 na*e 3e 


^ cos i ^ 1 ^ 

dt na* (1-e*) ^^*sin~i" 3u na* (1-6*7 sin i 3fl 


d» , _1 ^ , 

dt na* (1-e) '^* sin i 31 

M _ 1-e* iX ^ ^ ^ ^ 

dt na*e 3e na 

I- 

1/2 • 3 / 2 

where n is the mean motion (u a ' ). 
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The equations o£ motion can be approximately integrated 
under certain assumptions. Because o£ the smallness of the effects 
of tesseral harmonics, they can be treated (to first order) as 
linear perturbations about the orbit produced by only the central 
force and the secular second zonal harmonic term [£mpq *= 2010 ) 
in the potential. Under these assumptions, Kaula (1966, pp. 40, 49) 
gives the solutions foa* the ajiiplitudes of the perturbations of 


the Kepler elements 


due to each harmonic component V 


V.mpq 


as 


iimpq c ‘^[CA-2p)w+Cil-2p+ci)M+MCi'i-0)] 


Ae. -„a‘ ‘^n^-2p^g) - (t-2p) 1 


'ilmpq'^'^^e 


na"' -"e [ C Jt- 2 p) (. 1 - 2 p+q) e) ] 






Umpq 


■cot 


na* ‘"[CJt“2p)w+(£-2p+q}M+mCfi-e)] 


Ai 


„ i -..1 J,„ , 

^mpq e na^'^^Cl-e^) ^^^sin i[(il--p}.li+(3i-2p+q)M+inCiV-e)i (35 
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C ^ s in i [ Cii “ 2p) Ji+ ( Jt - 2p+q)M+m(C^- 6) ] 




^ [-Cl-eMe"H3Gj,_yae)+r,(An)G. IF. 
AM«„__=ya": ~,rr^ ^ ?m 


^mpq ®na^'*'-^[CJl-2p)i;)+()l-2p+ci3M+inCW-0)] 


[ (Jl- 2p) (i+ (a- 2p+q)M+m(fi- 0) ] * 


When examining these expressions, one can see that 
under certain conditions the frequencies C'l^) in their 
denominators can go to zero: 

^Umpq = C^-2p)w+Cii-2p+q)M+mCft-03 = 0 (4) 


This is known as the resonance condition. When this happens, 
one has exact commensurability between satellite motion 
and the earth's rotation yielding a perturbation from the 
longitudinal dependent terms of the geopotential analytically 
approaching infinity. This is known as deep resonance . 

Of course, other forces are acting on the satellite, so 
that the orbit usually simply passes through this condition 
of perfect conuaensurability , Typically, atmospheric drag 
is the dominant other force. 
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However, for a very large number of satellites, a 
situation exists where the resonance condition is only 
approximated, yielding substantial perturbations on the 
satellite nevertheless (Wagner and Douglas, 1969). This 
happenstance is called shallow resonance ^ These effects 
start becoming a problem for orbital operations and 
precise orbit determination wlien the resonant period, 
whidi is 

1 1 

TT ® ' 

tp (£“2p)w + (£- 2p+q3M+m y) 


starts to approach a few days duration. Unfortunately, in 
shallow resonance, all the frequencies of a given order 
satisfying Equation (4) are almost the same. It will be 
difficult therefore to distinguish these effects over a 
short time. The problem is to separate the information 
in order to determine specific gravitational harmonics. 
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ANALYSIS OF SHALLOW RESONANCE 


n " fi ® <iu *-M+m (Q- 0) , (m H ill revolutious/ 
jii ^ in ^ u 

day) be the dominant resonant longitude, and E be any o£ 
the Kepler elements. The effect associated with this 
longitude will be of order one compared to the fringe 
resonances, with longitudes q+w, which are 

of order lOe (Allan, 1973, p. 22"4) . Note that in one 
period of i|/ «, to will be essentially constant if 

• 9 lU ^ U ^ 

<|»m 0 usually is for shallow resonant orbits. 

As long as the resonant period is less than about 
100 days, Kaula ' s (1966) linear perturbations , (equations (3)) 
are valid (Gedeon, 1969) : 

^ ^0 X] '^Ampq ^^S '^'^mpq f 

Relevant ^ 

Umpq 

where 


*^Jlmpq ■ (5- lp)w+(£-2+q)M+m(«-6) , 
ignoring the phase as in equation (1) , 
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even 


for E ■ w,n,M 


L’^tm »"^tm 


i*m odd 


^^C,S “ ^^impq 


(7) 


^£m * ^£m 


l-m even 


for E ■ a ,e , I 


-■^£m ' ^ImJ 


£-m odd 


Note that (C,„,s,„) - {cos m sin m X^J and iE,„p^ 

(given by a right-hand side of (3) without is inversely 

* * 2 

proportional to 4'£„,pq (or for AM) which is small in 
resonance. 

For the five principal resonant frequencies q, 
*m,-l» ^m,^l» %,-2» '^m,2» (6) becomes simply: 


h * 


2 *n,0 • ♦m.fli 

I RESONANT, 
q“0 


[AEf, cos ^ 


£ RESONANT 

q»“l 


^ [AEj, cos ^'^m.ol^'**)] • 

£ RESONANT 

q»“2 
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In Equation (K) the resonant t's for the q*^l terms 
are always even and the resonant t's for the q«0 (dominant) 
and q"^2 terms are always odd. 

Expanding the cos and sin terms in (8) and collecting 
terms in cos sin * cos u and sin u, cos 2ut and 

m f u Ml g w 

sin 2<t>: 



♦ cos 2b) 


AEp * 2u/* AEq * 

t RES ^ t RES ^1 I IRES ^ t RES 

q «2 q --2 / \ q »2 q "*2 / 


♦ sini|/ 


m,0 


AE, 


t RES 
q-0 


* sinu) / ^ ^ 

I t RES 
\ q*l 


AE, 


* cosu/ AE- * 

It RES ^ I RES ^1 

\ q»l q«-l / 

t RES V 


cos2b) / AE--^ 

AEg\ 

^ sin 2b)/ ^ ^ AE- 


It RES ^ 

t RES 

\ t RES 

t RES I 

V q»2 

q»-2 ^ 

\ q«2 

q- * 2 ^ 
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( 9 ) 



Hquation shows that the (lumped) aoefficioiits (for each 
element) of the q and q terms (determinable in 

one period of q) are themselves sinusoidal functions of 
a slowly varying argument of .ierigee; 


C* C,, + G COS6J + C. sinto + G-„ cos2w + G,„ sin2w + ... 

(103 

S* = S, + S rosy + S„ siny + cos2w + sin2w + ... 

U ^ ^ ^ 4i^ 


The components of the lumped coefficients (C*, S*) depend 
on the resonant geopotential coefficients, as well as the 
a, e and i of the orbit. 


The relation between the lumped coefficients and the 
actual tracking infoi*mation is straightforward, but tedious 
to write out in detail. Essentially it is the same as the 
relation between the tracking information and the orbital 
elements. If 'O' is a tracking observation, it is clearly 
a function of the orbital elements E. These in turn are 
given (to first order) by; 


E “ Ey + E* 


cosij) 

siniii 


m,0 
in, 0 


over as many lumped coefficients E* as are necessary to 
describe the variation. This suggests a simple scheme for 
determining all the lumped coefficients for a single tracwing 
arc. Tracking residuals in 'O' can be resolved by differential 
correction to E^j and E* through the condition equations; 


r 


30 

COST' A 

m , 0 

AO = — 

AEq + AE* 

3E 
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The cbsei’Vation ’partials* yO/SH are laborious to 
calculate but are readily available in existing differential 
programs [Lerch et. al. 1974J, This is analogous to the 
method chosen by Riegber (19 73J to deterniine shallow resonant 
constraints. Riegber however collects a boundary (rather 
than an initial) value solution to the orbit. He also appears 
to make no use of the known frequencies of the problem, pre» 
ferring a general Fourier analysis of the resonant elements. 

However, the principal periods in Reigber's solutions are the 

« 

full arc lengths, chosen as q. His GEOS 1 1 solution for 

these are compatible with our results, as will be shown later. 

One of the goals of our analysis, however, is to identify 
the dominant information in current fields. Where possible, we 
want to uSi vHe good geopotential solutions, already developed 
at great tApeiise, Our central liypothesis is that the resonance 
information is almost entirely in the along track (aw+aM+aScosi) 
variation, To give a concrete example, consider the orbit of 
GEOS II (a » 1,209 e.r., e .033, i = 105, 8^^}. The root sum 
of squares of all resonant terms contributing to this orbital 
component is indeed 95% of the total perturbation. The 
information in the semi-major axis variation is almost the 
same as this since its integral controls the along track change 
in resonance. Table 1 give the resonant ’a’ perturbations on 
GEOS II from x 10‘^/£^ to (30,13). 


Clearly, the q=0 terms dominate but the q=+^l terms are 
also significant. To illustrate the development of the lumped 
coefficients, the quantities A (a) are the components or 

AEg (for E=a) for coefficients or of /? x 10’ in 

equation (9), where AE^ g(q=+l) = ^£m^ 

= ’ Therefore, to find the partial contribution 

(or sensitivity) due to each (unknown) coefficient (AE) the 

perturbations in the table must be divided by /I x 10 /£ . 


11 


In fact they can also be adjusted by a constant (as a set) 

without changing the relative information content of the 

terms in the lumped coefficients C* and S*. Multiplying 

2 

each term of Table 1 by A /400 gives the partial contri- 
butions presented in Table 2, 

To obtain non-dimensional sensitivities, these partials 
are divided by the maximum contribution (from t=15, q»0) 
yielding the sensitivities presented in Table 3. 

Thus, from this final sensitivity table the lumped 
coefficients determinable from observation of just the *a’ 
variation are (to q=i2 terms) 


(Cp,Sg) = -.886(C,S)^3^j. + 1.000(C,S)^5^j^3 


+ .456(C,S) - 

.02orc.s)ig_i3 - 

*1®1(C»S)2x,13 

•.156(C,S)23^^3 - 

.076(C, S3 25^13 - 

0 . 010 (C, S) 27 X3 

+.022(C,S)2g^j3 + 

» • • 

(UJ 

-.514(-S, 0^4^13 . 

.034C-S.C)j^^j3 


+ . 160 ( -S , L) ^3 + 

.U0C-S.C32o_i3 


+ . 0 3 5 ( - S , G ) 2 2 ^ 3 

.011(-S,C)24^23 


-.025(-S,C)25^13 

.021C-S,C)2g^TL3 


- .011(-S,C)3Q^^3 . 

• « 

(12) 
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^5*®S^ " * ^ 14 j 13 ^ *^^^^^'^^16,13 


+.041(0,5)28,23 ' -068(0,5)20,13 ' 

.092rC,S)32^33 

-.064(0,5)24,23 - .022(0,5)20,23 - 

•°“SCC.S)28_33 


+ .021(C,SJ2q^jI^2 + • • • 

(Gc2>®c 2^ “ ” * ® ^ 13 ^ 13 ^ 15 j 13 

+.012(^,5)27^23 ^ ^ 19 , 13 

+.006(0,5)22 23 ^ .005(0,5)22 13 

+.002(0,5)25^23 ' * »S) 27 ^23 
“.003(0,5)20 13 ^ *'* 


( 13 ) 


(14) 


^^52 »^S 2 ^ 


+.031(5,-0)22 23 ^ 
“. 023 ( 5 ,- 0 ) 27^23 
“. 015 ( 5 ,- 0 ) 22,23 " 
+. 007 ( 5 ,- 0 ) 25,23 " 

+ . 006(5,-0)29,23 ■*■ 


. 010 ( 5 ,- 0 ) 25,13 

. 027 CS,-C)„_i 3 

. 001 ( 5 ,- 0 ) 22,23 

. 008 ( 5 ,- 0 ) 27,23 


(15) 


One can also repeat the above computations using the 
sum of the perturbations Am, AM and cos i•A^^ to derive 
influence coefficients for observation of an effective 
along track variation: 
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CCo.Sq3 “ 

872(“S,C3j^3^j^2 ■*■ 

l-0»(-S,C)is,i3 * 


-.020(-S,C3iy^l3 - 

.189(-S,C)2i_13 - 


-.U82(-S,G325^j3 - 

, 011(-S ,C3 27 ^13 ■*■ 

(Cc.Se) “ 

-.516C-C,-S3^4^j^2 + 



+ .11G(-C,-S32q^]^3 + 

.039(-C.-S322,i3 


-.027(-C,-S32o,13 “ 

,024(-C,-S32g^j_3 

(Cs,Sg3 = 

.17K-S, €3^4^13 

.189(-S,G3i,^13 " 


-. 073(-S,C)2q^2.3 * 

.098C.S,G322,13 ’ 


-.025C-S,J326,13 " 

. 01D(-S, 0328 ^13 

(Cj 22 »S(223 

= -.040(-S,C)3^3^^3 

+ .020C-S, 03^5^13 


+ .014 C"S»C3^7^^3 

+ . 008C-S, 033 ^ 0 ^ j ^3 


+ , 00 7 (- S , C3 2 ]^3 

+ »006(-S, 0323^13 


+ .002(-S,C3 25 ^;L3 

- .OOK-S, €3 27 ^15 


-.004(-S,C32ci^13. 

• • 

(*^52 >^S 2 ^ 

= + , 033 (C ,S3 13 13'*' 

.009(C,S3i5^13 


”.024(C»S3^^ j ^3 * 

• *^“^^^'^^19,13 


-.016CG, 8 ) 21^23 

•001(^.5)23^13 


+ . 007 CC,S3 25 ^j ^3 + 

.009(G, 5)27^13 


(1G3 


C17] 


(183 


(193 


+ . 007(CjS32g 


(203 
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TABLE 1 

PCBTURBATION AMPLITUDE IN • ^U) matwi 
(m - 131 





q^JL 

qr *1 


13 

-.011 


-2589 


-.103 

14 


-.448 


-.847 


16 

034 


2.196 


.011 

16 


.214 


-.149 


17 

-009 


.779 


.030 

18 


.1M 


.091 


16 

-.014 


-.026 


.023 

20 


.026 


.109 


21 

006 


-203 


.012 

22 


- 029 


066 


23 

-002 


-.146 


.003 

24 


-.03? 


023 


26 

004 


-.060 


-.002 

26 


-.017 


-.001 


27 

003 


-.007 


-.003 

28 


- 004 


-.009 


29 

001 


013 


-.003 

30 


003 


-.009 


RSS 

.04 

0.61 

3.2 

0 86 

.11 




TABLE 2 





PARTIAL CONTRIBUTION ^ la) 

m 





m - 13 



J_^ 

q‘ -? 

-1 

q ■ 0 


9 " *2 

13 

-.006 


-1 094 


- 043 

14 


-.220 


-416 


16 

.019 


1 235 


006 

16 


.137 


-.096 


17 

-007 


.563 


.022 

18 


.123 


.073 


19 

-.013 


-.024 


.021 

20 


.026 


.109 


21 

-.006 


- 224 


.013 

22 


-.035 


079 


23 

002 


-.193 


.004 

24 


- 046 


033 


25 

006 


- 094 


-.003 

26 


-.029 


-.002 


27 

006 


-.013 


-.007 

28 


-.006 


-.019 


29 

002 


028 


-.006 

30 


-.006 


-.020 
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TABLE 3 

SENSITIVITY COEFFICIENTS 
FOR GEOS-II RESONANCE 
^E(a) 




-J 

13 

-.004 


14 


-.178 

16 

.016 


16 


.111 

17 

-.006 


18 


.100 

19 

-.010 


20 


.021 

21 

-.004 


22 


-.029 

23 

.002 


24 


-.037 

26 

006 


26 


-.024 

27 

004 


28 


-.006 

29 

.001 


30 


♦ 005 


a- 0 

Q ■ »1 


-886 

-.336 

-.036 

1 .000 

-.077 

.006 

.466 

069 

.018 

-.020 

.089 

.017 

-.181 

.064 

.011 

-.156 

.026 

.003 

-.076 

-.002 

-.002 

-.010 

-.016 

-005 

.022 

-.016 

-.006 
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As predicted, these sensitivities only vary slightly 
from the values derived for Aa, 

Similarly, the information content from the element 
variations in ’e ' and 'i' is almost the same or predictable 
from the 'a* variation since: 


rCl-q) cos i-m] 

A i/ Aa ~ — TT7 f and 

2a(l-e^) sin i 


2 1/2 2 1/2 

Ae (l-e^3 [(1-e^) -(1-q)] 


Aa 2ae 


( 21 ) 


Note that these ratios depend on the frequency q as well as 
the resonant order ra, but they are independent of degree. 
Therefore, the sensitivity tables for e and i information 
are the same as those for 'a' with the following adjust ~ 
ments (Table 4) : 


TABLE 4 

FACTORS OF THE 4» TABLE FOR 
GEOS II OBSERVATIONS /le AND 



q « -1 


-13.54 


-1. 

AV= 4I/-13.54 

1. 


q** 0 1 

-13.27 -13.00 

-5.45 X 10”^ +1. 

.960 .960 


18 



Table 4 shows that the sensitivities for Ai observations are 
virtually unchanged from those for ^a while the Le information 
appears to be significantly altered. But in fact, the Ae 
information is almost entirclv in the q«+l terms which are 
predictable from the 'a' variation, In fact; 


C^Ue) = SgC^a} 

SgCAe) = C^Ua) 

(22‘J 

CgCAe) - -S^CAa) 

S^id.e) = -C^CAa) 

so that no new information is added by observations of 
eccentricity variation. 

In summary, we have hypothesised that nearly all the 
information in shallow resonance is contained in observations 
of either the semi-major axis or along track variation. As a 
consequence, analysis of tracking data within each short period 
q) may be made in terms of a simple "lumped” coefficient 
set. Such a set will vary sinusoidally with the long period 
of the argument of perigee, A similar analysis can be made 
for the higher order resonances of 2m, 3m, ... , on the 
same orbit. 

We have chosen the GEOS- II (1968 2A) satellite to test 
our method of analysis. GEOS- II has a principal resonance 
period which is approximately 6.5 days. This satellite was 
selected for three major reasons: 

1) It has been heavily tracked with very accurate 
instrumentation. 
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2) This satellite afforded us the opportunity to 
; calibrate our new technique for identifying 

resonant constraints since it was used in 
almost all recent global geopotential 
solutions , 

3) The satellite is largely unaffected by 
atmospheric drag. 

The total resonance perturbation on GEOS-II along 
track is approximately bOUm, and is an enormous effect 
when compared to the Im laser ranging and ~1VS0 camera 
instrumentational accuracy which was employed to track 
this satellite. 



ANALYSIS OF GE0S“1I DATA FOR RBSONANCE DETBRMINATION 


Following the hypothesis developed in the previous 
section, the entire laser and camera data set available on 
GEOS-II (ill IDoS-oS.) was divided into forty-one o.S day 
segments. The GEMl gravity model (Smith, Lerch and Wagner, 
1973) was then employed and the orbital state was estimated 
using this tracking data. We, however, removed all m®13 
terms from the GEMl model and recovered a lumped value of 
C(13,13) and S(13,13). These recovered values are plotted 
(fully normalised) against m in Figure 1, The results for 
the various GEODYN solutions are summarized in Table 6. 

The GEODYN orbit determination system (T. Martin, 
1972) was used for the resonance determination. GEODYN is 
a Bayesian least-squares, multiarc, multiple satellite 
orbit and geode tc parameter estimation system based upon 
Cowell type numerical integration techniques. Modeled 
parameters normally include luni-solar gravitational per- 
turbations, solar radiation pressure, atmospheric drag, 

EIH polar motion and UTl data and the GEMl geopotential 
model. Initially drag was modeled but not adjusted in 
the single arc solutions. Considerable correlation was 
found between the drag coefficient and the recovered reso- 
nant coefficients. Later, small but significant adjustments 
of drag were made in a multi- arc solution, having low 
correlation with the lumped coefficients. 


A W'eighted least-squares solution was then performed 
to find the values for the sinusoidal terms in equation (10) 
from the lumped values for (13,13) recovered from the 
GEODYN orbital analysis over each 6.5 day arc. This 
solution found the following values describing the lumped 
coefficients : 
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7b 13 

342 41 

0 1b 

018 

088 

12 

-90 24 

79 bO 

331 83 

0 48 

044 

080 

13 

-873b 

79 80 

323 41 

0 18 

022 

1 00 

14 

-92 87 

74 08 

312 42 

027 

0.31 

086 

1b 

98 07 

71.3b 

298 83 

034 

034 

1 28 

IS 

38 73 

7280 

280 07 

023 

022 

0 81 

17 

33 bb 

08 7b 

278.73 

0 21 

014 

087 

18 

-21 49 

96 80 

286 66 

0.23 

017 

1 00 

10 

-1884 

90 82 

26641 

026 

0 16 

1 11 

20 

-17 82 

43.31 

143 66 

036 

02b 

041 

21 

23 81 

4b 01 

134 27 

1 26 

084 

088 

22 

-20 80 

44 4b 

12b 13 

0 22 

023 

090 

23 

34 08 

48 88 

112 70 

066 

048 

086 

24 

-41 92 

48 78 

100 87 

0 28 

028 

084 

2b 

48 2s 

93.42 

81 48 

0 24 

033 

1 06 

2S 

-b1 84 

01 28 

81 96 

027 

0.37 

088 

27 

-b9 78 

58 33 

72 14 

0 62 

1 08 

085 

28 

-90 25 

94 44 

81 66 

0 40 

046 

09b 

29 

-70 10 

73.61 

62 40 

040 

036 

1 14 

30 

-18 18 

52 44 

246 74 

0 28 

016 

1 14 

31 

-13.1b 

48 80 

234 82 

0 13 

0 11 

1 64 

32 

-1040 

47 02 

226 68 

036 

0 16 

1 38 

33 

-1784 

48 87 

169 36 

0 18 

031 

080 

34 

-23.11 

48 42 

148 70 

0 16 

0.22 

1 01 

3b 

28 20 

54 88 

137 01 

0 18 

022 

1 34 

30 

-33 0b 

63 32 

128 82 

0 16 

0 24 

086 

37 

-1081 

62 78 

178 18 

OM 

026 

1 16 

.'8 

45 08 

59 17 

111 52 

0 18 

013 

1 01 

39 

-41.70 

b1 83 

101 08 

033 

023 

1 28 

40 

58 17 

66 28 

81 74 

062 

0.37 

0'7 

41 

-89 80 

86 61 

81 88 

066 

068 

0 77 


23 


original paob b 

OF POOR QUAUni 



Gq « -40.UU86 X 
Sy " 61.2619 X lO"^** 

= 15.6332 X 10“^ 

C*. =» -29.0172 X 10”^^ 
c' « -11.3765 X 10'^ 
s[ = - 2.8247 X lU"*^ 


Since no other resonant term is modeled the (C,S )^2 2^2 
values theiiisel\?es can be considered the lumped coefficients in 

t I 

in equation (10). We have labeled these C ,S . In equations 
(11- 20) , the full model for these coefficients are (arbitrarily) 
normalized with respect to the (15,13) term which has the 
greatest influence. 

The assumption here is that the resonant information 
is entirely along track (equations 16-20), Therefore, 
according to the convention in these equations (C*.,ev') 

I t 

- .872 (S ,-C ), The terms of these lumped resonant co- 
efficients are thus : 


^0 

= 

34.922 

X 

10-^ 

"‘"0 

= 

53.475 

X 

io-« 

s , 


25.328 

X 

10-^ 

c 

-C 

c 

= 

13.646 

X 

lu-s 


= 

9.930 

X 

10'^ 


s 

-2.466 

X 

io"“ 
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These coefficients can then be used in equations (16-18) 
to produce any three resonant coefficient sets modeling the 
three distinct frequencies for the q«-l,0 and +1 terms. Such 
a set is presented in Table 7, The q''+_2 frequencies were 
barely detectable in the lumped coefficients and were not 
successfully recovered from this data. 

This same analysis can be performed using the ^a 
constraints found in equations (11), (12) and (13), The 
same resonant coefficient set for Aa are also presented 
in Table 7. Since the two sets of constraint equations 
are nearly the same, these sets are also. 


TABLE 7 

RESONANT COEFFICIENTS 
FROM GLOBAL * GEOS-H TRACKING DATA 
(SINGLE ARC ANALYSIS) 



from 

from 


(d Along Track) 

(/I a ) 

C(1S,13) 

34,922 X 10~® 

35.438 X 10~® 

S(1S.13) 

-53.473 X 10*~® 

-54.263 X 10“® 

C(14,13) 

27,976 X 10~^ 

28.880 X 10“® 

S(14,13) 

47,408 X 10“® 

48,204 X 10”® 

C(16.13) 

27.228 X 10"^ 

29.330 X 10~® 

S(16,13) 

-29.845 X 10-® 

-27,204 X 10”® 


*Data ever a full rotation of perigee. 
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ERROR ANALYSIS 


The ORAN program CC.F, Martin, 1970) was used to 
perform a comprehensive error analysis of our resonant 
harmonic determination. In particular we wanted to know 
why the single arc lumped harmonics in 1969 were sys- 
tematically displaced from those in 1968 (see Figure 1). 
The ORAN program calculates the effect of unsolved- for 
(and poorly determined) parameters on the resoncnt 
determination. ORAN does this (without lengthy simula- 
tions) by computing numerical measurement partials with 
respect to a large number of unrecovered effects. Three 
kinds of problems were investigated: 

1. ORAN was used to perform a classical error analysis 
and gave both accuracy assessments and located the 
dominant error sources affecting the recovered 
resonance terms. The modeled errors included; 

• atmospheric drag at 40% error in a ballistic 
coefficient, Cj^: 

• tesseral and zonal harmonic errors at 25% of 
the difference between two independent gravity 
models (Martin and Roy, 1970); 

• tracking station coordinate uncertainties in 
the adopted set of positions employed for 
these solutions (Marsh et. al., 1973): 

• 


an error in u of 1 ppm. 


The results o£ this error analysis are summarised 
in Table 8, This analysis indicated that atmos- 
pheric drag was a large error source and could 
substantiallv bias the resonance recovery from 

r * 

eacli arc. This was surprising since the overall 
effect of drag on the rather high GEOS II orbit 
is small. Greater effects were seen, however, 
in the high solar cycle years of 1968-69. In 
fact, a large part of the discrepancy between 
the 1968-1969 data in Figure 1 could be due to 
drag error. Tlie non-resonant geopotential error 
is probably overestimated in Table 8 since the 
field used for the GEOS arcs (GHM-1) was deter- 
mined with mucli of the same optical data. Also, the 
gravity model error magnitude was scaled vo the SAO 
Standard Earth II Gravity Model (Gaposchkin and Lambeck, 
1970) whidi has been shown by Klosko and Krabill (1974) 
to yield approximately twice the orbital errors 
as that of GEM-1 on GEOS II orbits. Attempts 
were made to determine a drag coefficient, Cj^, 
in each of tlie 6.5 day orbital solutions, but 
excessively high correlations between Cj^ and 
the resonance terms (at times exceeding .90) 
prevented our having much confidence in the results. 

The ORAN program A\?as also used * v^erify the 
analytical development and po^ ^tion of the 
previous sections of this report. Values for 
all of the 13 (through degree 21) order reson- 
ance terms were modeled as error sources at the 
magnitudes listed in Table 9 as a and thereby 
liad their perturbations propagate into the re- 
covered values for (13,13) , ^ at the magnitude 
listed as B, By comparing quantities a and 8 we 
were able to numerically determine values fur 
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TABLE 8. EST,.* ATED ERROR IN RECOVERED (13,131c S 
VALUES FROM ERROR ANALYSIS FOR TWO SAMPLE GEOS-II ARCS 
(NORMALIZED VALUE X ;0 *) 


ERROR SOURCES 
(Magnitude) 




u 

(1 »>pni) 

DRAG 
140% OF 

Co> 

GRAVITY 
MODEL 
ERROR 
(.25 APL- 
SmO M-1) 

STATION 
COORDINATES 
(AT MARSH 
ESTIMATED 
UNCERTAINTY) 


ARC 

EPOCH 






RSS 

TOTAL 

680307 

C(13.13) 

.015 

3.068 

0.771 

a33^ 

3.181 


S(13.13) 

-.011 

0.464 

1.238 

0.929 

1 616 

690207 

C(13,13) 

.027 

-1 660 

3.073 

a747 

3.572 


S(13,13) 

.029 

2.579 

-0.415 

0.265 

2625 



TABLE 9 

ORAN CALIBRATION OF GEOS-II SHALLOW RESONANCE STUDY 


CO, SO 

lai VALUE 
PROPAGATED 
INORMALIZEO) 
X 10-8 

1^1 VALUE 
OF COEFF 
ABSORBED 
BY 113,131 
X 10“® 

ORAN 

ESTIMATED 

SENSITIVITY 

1 ^} 

— SENSITIVITY VALUE PREDICTED FROM THEORY 
SEMI FOR q- - 1.0.*1 TERMS 

MAJOR ASCEND. ALONG 

AXIS NODE ECC. TRACK 

CI15.13) 

4.444 

-5 174 

-1.164 

-1 129 

1 564 

-1.129 

-1.146 

SI1S.13) 

4 444 

-5.144 

-1.158 





C(17.13l 

3 460 

-1.840 

- 532 





S( 17,131 

3 460 

-1 842 

- 532 

— .515 

-1 182 

— .515 

- 530 

Cl 19. 131 

2.770 

♦ .087 

031 





SI19.13I 

^770 

+ 076 

027 

.022 

—1.726 

.022 

.023 

CI21.13I 

2.628 

* .514 

227 





SI21.13) 

2628 

♦ 507 

223 

.204 

-0 954 

206 

.217 

ICC, SCI • cot u 








-S(14,13) 

SI. 020 

4 160 

292 





CI14.13I 

SI 020 

4 106 

288 

.162 

- .000 

-102.83 

.165 

-S1 16.131 

39 063 

♦ 0 292 

- .027 





CI16.13) 

39 063 

-0.195 

- .018 

- .011 

266 

-122.50 

- 009 

-SI18.13) 

3 086 

♦ 0 086 

- 099 





CI18.13I 

3 086 

-0081 

- 094 

- .C50 

133 

— 26.71 

- 052 

-SI20.13) 

2.500 

♦ 0 046 

- 066 





CI20.13) 

2500 

-0 047 

- 064 

- .C35 

- 032 

43.96 

- 037 

ICS, SSI - tin CJ 








CI14.13I 

51 020 

-7636 

- .156 





SI14.13I 

51 020 

-6454 

- .132 

- .134 

.233 

1150.3 

- .188 

Cl 16.131 

3" i ' H53 

-7 539 

- 201 

- 204 

- .162 

- 75 32 

- 208 

SI16.13) 

39 063 

-7530 

- 201 





CI18.13) 

3 086 

-.1373 

-.0463 





SI18.13) 

3 086 

-.1597 

- 0539 

- 044 

— .605 

-357.2 

- 042 

CI20.13) 

2500 

♦ .1756 

♦ .0731 





SI20.13) 

2 500 

♦ .1600 

♦ 0667 

.073 

- .542 

244 9 

080 



the sensitivity coefficients, cliecking tlie analytic 
results. Table 9 presents these results. As antici* 
pated, the along track constraint best represents 
the relative sensitivity of the 13^^^ order terms for 
the q-0, and composite q=+l frequencies for GEOS- II. 

The fringe terms are less well modeled than t)ie dominant 
constant and terms, especially the small cosw 
terms for this sa.nple u « 7 5" orbit, 

3. ORAN simulations were also performed in a manner 

similar to the preceding, but simulating the ad- 
justment of more than one pair of 13^^ order terms. 

This analysis was done to assess how much resonance 
information was contained within each b.S day data 
set. By propagating the effects of the unsolved for 
and therefore neglected resonance terms into the 
orbit, and getting an estimated satellite positional 
error, we found tliat when two even and two odd degree 
IDairs of coefficients were recovered, the remaining 
orbital ei’ror from all other resonance terms was 
estimated to be les ■ than ■ j cm. When two even and 
a single odd pair of coefficients are recovered, 
the estimated orbital errors were at timfs, 2 meters. 
With a single odd pa^r recov^ered, the estimated orbital 
errors due to neglected resonance, at times were 10 's 
of meters. We therefore deduced that two even and 
two odd pairs of coefficients could recover the total 
resonance information for GEOS-II. 


The problem of atmospheric drag errors was still present and 
we therefore decided to recover a single set of two odd and 
two even pairs of resonance coefficients in a multiple-arc 
solution using 14 6.5 day orbits well distributed over the 
apsidal period. Two pairs of 26^ order day period) 

terms were also estimated. In each of these arcs, a was 
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independently estimated. However, since all 14 arcs contributed 
information to the resonance recovery, the correlation between 
the recovered and the resonance terms was satisfactorily 

reduced to no more than 0,0, and seldom exceeded 0.3, 

It is these results which we have adopted as best for 
t h. 

this report. These 13 order coefficients are presented in 
Table 10 and are used for the theoretical resonance values 
plotted in Figure 1. 

C o m oarisons with Compreliensive Gravity Model 

Solutions 

There are many comprehensive grav^ity models which 
liave been produced using satellite tracking data. Some 
used data from GEOS-Il, while others did not. One can 
get some estimate of the consistency bet\ceen various models 
and also compare tlie results inferred from these models 
for GEOS -II with the results w»e have obtained using our 
numerical analytical technique. 

By taking equations (10) through (20) derived from 
(dfjj •<- uM + cos i'Afi) and substituting coefficient values froiu 
a given gravity model, one can compute a value of the lumped 
coefficients for eacli of the gravity models. Table 11 
presents these results. 

Table II indicates very good agreement between our 
GEOS- II multi-arc analysis and the results obtained from 
comprehensive geopotential solutions wdiich had a strong 
presence of GEOS- II data. The lumped coefficients for 
the comprehensive models are remarkedly consistent (on 
the whole) in spite of a fairly wide divergence of actual 
coefficients (see Figure 2 and Table 12). In Table 12 we 
have listed b satellite only solutions (and their GEOS- 


TABLE 10. QEOS-II RESONANCE HARMONIC RECOVERY 
USING A MULTI-ARC SOLUTION 


ALONG TRACK 

COEFFICIENT NORMALIZED CONSTRAINT LUMPED 

VALUE H 10*® VALUES 


C(13.13l 

-65.166 



S(13,13) 

70 009 

- 

38 1376 

C(14,13) 

27.064 

Co- 

>^55 9037 

SI14,13I 

53.582 

Cc - 

♦13 3597 

C(15,13) 

-18.744 

®C " 

♦27.9689 

S(16.13l 

5 206 

Ss‘ 

8 5028 

C(16,13l 

20 516 

Cs- 

-6 8818 

S(16.13) 

-11.993 

■ 

U 

2.7169 

C(26.26l 

-63.160 

®2C“ 

2 2488 

S(26,26) 

23.582 

C2S- 

- 2 3294 

C(27.26) 

- 3.035 

**2S“ 

2.3611 

S(27,26) 

30 581 
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TABLE 11 GRAVITY MODEL COMPARISON FOR RECOVERED 
LUMPED GEOS-II REIONANCE COEFFICIENTS 
(NORMAL ZED X 10®l 


LUMPED 

COEF 

(M 10*1 


GEOS-II 

MAGNER 

KLOSKO 

ANALYSIS 

(MULTI-ARCI 


GEOS II 
WAGNER 
KLOSKO 
INDIVIDUAL 
ARCS 


MODELS NOT USING GEOS-II 




MODELS liSING GEOS-II 



RAPP ytONOULIS 

DOUGLAS* 

SAO 

SAO 

GEM1 GEM4 GE1AS 

GEM6 

PCS 30PPIER 


50 II III 62 





(19671 

119731 

119681 

(19691 

(19701 

(19741 

(19711 

(19721 

(19741 

(19741 

(19741 


*o 

38 138 

34 92 

17.41 

21 67 

18 10 

18 10 • 

37 78 

29 96 

37 52 

40 06 

39 93 

40 03 

38 71 

3*88 

Co 

SS90« 

53 47 

21 41 

6287 

68 04 

68 04* 

56 57 

62 48 

56 10 

55 73 

55 02 

55 02 

55 32 

54 54 

*c 

27 989 

CoiUl 

25 33 

1482 

2-32 


48 18 

4 77 

27 54 

29 87 

27 48 

28 33 

28 16 

27 22 

27 07 

Cc 

13 360 

1365 

10 37 

4 02 


4 22 

1806 

16 29 

11 32 

13 59 

12 97 

12 94 

15 22 

16 64 

H 

8603 

Stnui 

993 

3 39 

93S 


1 38 

11 43 

701 

964 

5.59 

598 

606 

7.77 

824 

Cs 

-6 882 

-247 

-465 

- 2 95 


-15 77 

- 789 

- 5.41 

- 526 

- 2.43 

- 3.53 

-350 

-3 74 

-352 


2 717 

Cot 2UI 

— 

— 

— 

— 

— 

357 

3 43 

306 

281 

2 74 

280 

2 49 

306 

s 

2 249 

— 

— 

— 

— 

— 

784 

2 37 

1 65 

1 28 

1 70 

166 

1 79 

200 

s 

2 329 

Swi 2U 

— 

— 

— 

— 

— 

-4 47 

-335 

-1 20 

- 1 18 

- 1 31 

-1 35 

-1 72 

- 3 01 

s 

2 361 



— 

— 

— 

— 

456 

3 97 

334 

1 80 

2 49 

255 

213 

225 


•used YIONOULISIISeSI OOO DEGREE TERMS 
*FROM ALONG TRACK CONSTRAINTS 




TABLE 12 13 

’'^.TBDER SATELLITE SOLUTIONS 

UNITS. 10 

-8 




GEM 1 

GEM 3 

8CS 62 

UOA9' ER 

RIEGBER 

GEM6 

10 V* 

AVQ 

RMS 

IX 1 

flML— 











10^/7^ 

IX 

«3 3 

-26 2 

698 

82 

444 

68 8 

68 2 

62.8 

1X3 

226 

m 

70 4 

96 1 

68 3 

70 

68 1 

68 9 

68 2 

71 8 

11 2 

19 0 

14C 

174 

X2 

24 8 

32 

36 3 

14 6 

61 0 

26 7 

78 

14 8 

US 

4«4 

46 

42 0 

48 

466 

264 

61 0 

368 

16 7 

MS 

IK 

-2S4 

- 1 4 

-26 1 

-31 

-13 1 

-22 6 

44 4 

-19 9 

99 

22 3 

Ibt 

- 06 

12 0 

- 26 

-11 

- 84 

- 22 

44 4 

- 1 9 

74 

16 6 

IK 

IS 3 

4 4 

66 

18 

37 4 

09 

38 1 

138 

123 

31 4 

16S 

- 34 

-21 9 

- 26 

- 4 

-178 

-13.2 

38 1 

10 6 

76 

19 6 

17C 

6S 

32 4 

13 1 

19 

10 4 

104 

346 

16 1 

8 7 

26 1 

17$ 

111 

43 6 

23 8 

27 

190 

18 8 

346 

239 

10 1 

282 

ISC 

- 8 7 

- 39 

-203 

-10 

64 

-22.6 

M8 




1SS 

38 0 

89 7 

-39 8 

-29 

-20 1 

-61 9 

M9 




IK 

-28 6 

- 6 2 

-16,3 

6 

B2 6 

-244 

27 7 




I9S 

464 

- 03 

-23 3 

-:i 

-200 

-300 

27 7 




TOC 

2 1 

27 7 

6 9 

19 

366 

- 20 

260 




70S 

1 6 

-72.6 

1 1 

14 

2 1 

-31 7 

26 0 




21C 

-280 

-188 

-26 3 

- 9 

466 

-27 0 

22 7 




71S 

1 0 

26 3 

16 2 

- 1 

19 2 

10 8 

22 7 




2X 

41 3 

-ISO 

-226 

• 21 


41 2 

18 9 




22S 

60 

39 8 

16 7 

20 


-169 

18 9 




23C 



-11 4 

-19 







23S 



24 

-17 







24C 



-22 8 

10 







24S 



- 4 2 

- 1 







2K 



11 6 

-13 







26S 



78 

-19 







2K 



8 6 

11 







26S 



- 06 

- 4 







27C 



- 3 6 

-16 







27S 



128 

2 







28C 



66 1 








28S 



- 8 7 








2K 



-164 








T9S 



-16 9 








30C 



200 








30S 
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Co 

66 10 

66 76 

66 32 

64 64 

63 61 

66 02 

70 2 

64 9 

71 

1 0 

*o 

37 62 

40 09 

38 71 

37 88 

40 48 

39 93 

70 2 

39 1 

1 13 

1 6 

Cc 

11 32 

1267 

16 22 

16 64 

11 86 

12.97 

27 0 

136 

1 86 

6 9 

»c 

20 87 

27 46 

27 22 

27 07 

2767 

28 33 

27 0 

27 9 

96 

36 

Cs 

- 5 26 

- 2 60 

- 3 74 

- 3 62 

- 3 71 

- 3 63 

11 9 

- 3 7 

81 

68 

s$ 

964 

6 69 

7 70 

8 24 

1067 

6 98 

11 9 

80 

1 82 

16 3 
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NORMALIZED COEFFICIENT VALUE X 10' 


FIGUPE 2. 

COMPARISON OF GLOBAL GEOPOTENTIAL SOLUTIONS 
FOR 13th ORDER TESSERAL HARMONICS 
COMPAREO TO FOURIER TERMS 
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lumped coefl'icients) with strong GEOS-II tracking. Gllfll (.bmith 
ot.al, 1973] contains optical data only, GEM3 ILorch et.al., 

1972] has electronic (Tranet Doppler, C and S-baud radar and 
laser) and additional optical data added to the GEMl tracking. 

The data is employed at full weight (according to the accuracies 
o£ the systems as judged by the arc residuals). Thus, GE0S>I1 
tracking dominates the GEM3 geopotential, but this Jteavy track- 
ing had some deficiencies which were remedied in GEMS and 
later solutions. The chief deficiency was the poorly known 
Doppler stations. In GEMS [Lerch et.al., 1974] the electronic 
and optical data for the GEOS (I and II) orbits were down- 
weighted to reduce the effect of these processing errors. 

With PGS62 (F.J. Lerch, Private communication, 1974), the 
station positions and data biases were resolved and the 
GEOS -II data (with additional Doppler and laser tracking) 
had full weight again. 

In Table 12 the Doppler solution uses only Doppler 
data on 9 distinct satellite orbits with heavy GEOS-II 
coverage. The lUegber and Ilk (1975, Table 2, Col. 3) 
solution uses optical and laser data reduced to 5b (resonance) 
condition equations on 6 orbits including GEOS-II . 

In spite of these differences in observation and 

data reduction, tae GEOS-II lumped harmonics are relatively 

stable compared to the geopotential itself. [In Table 12, 

we compare the RMS of coefficient variation with Kaula's 

rule (ll)''^/£“) . This is a uniform measure of precision 

over all degrees. For the lumped coefficients tlie estimate 
“5 

(10 /i") is actually the root sum of squares of all terms 

in each lumped component.] Except for S^, the lumped coefficients 
have a precision about an order of magnitude greater than the 
geopotential harmonics. The greater scatter of S values may 
be due (subtley) to the scarcity of GEOS-II observations 
around u=190®. This comparison is presented graphically 
in Figure 2. In general the lumped coefficients appear to 
represent excellent constraints for the GHOS-II tr? “?ng record. 
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Compurison with Dout*las» Marsh and Williamson 


Douglas, Marsh aii4 Williamson (19(39) porformod an 

analysis to recover resonance terms from timing errors 

(along track errors) in the Rosman G1^\RR data over a 5 .day 

th 

arc of GEOS-II. They modeled the odd degree, 13 order 
terms with the values of Yionoulis (190B) and attempted to 
enliunce this set of coefficients by solving for a value of 
tl4,13} Wfiich would combine to properly model the GfOS-II 
resonance perturbations. They assumed all the timing error 
was taken up in tlje resonance perturbation of the mean 
anomaly, Tliis is the second term of the expression for 
M in equation (3) . A more complete expression for the 
along track error [Aw+^m+cosi • which we use, differs 
from tlieirs by about lUO. In addition the Douglas-Yionoulis 
values apply only to the lij.ited argument of perigee during 
tiieir 5-day arc (~3SU"j . Nevertheless a comparison of our 
lumped coefficient with theirs at this perigee for GEOS-II 
is reasonably close. 

Using the constraint for along track information only 
[Table 11 in (10) for our (multi-arc) solution; equations 
(10-15) in llO) for tlie Douglas-Yionoulis set]; this comparison 
is presented in Table 13. The o 'global' constraints them- 
selves, computed from the Douglas-Yionoulis field are much 
poorer. But GEOS -II data is only represented in that field 
by a single satellite arc (see Table !)• 


Coiuparlson with Riegber's (1973) Constraint 

In his 1973 paper Riegber has shown how it is possible 
to derive resonance (or other periodic) variations from a 
'Tourier' solution of the satellite's motion as a boundary 
v:flue p>roblem. Instead of the 'natural' frequencies (?(() of 
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TABLE 13 

COMPARISON OF DOU6LA8-YIONOUU8 AND GEOS-li 
GLOBAL RESONANCE ANALYSIS FOR 
LUMPED COEFFICIENTS AT 
tJ - 349°.e6 



C X 10* 

S X 10* 

DchjoIii, ttal, 1969 

76.02 

66.26 

GEOS'-! 1 Rttcnanse 

Antlyiii (this rtport! 

67.81 

64.14 
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tile orbital perturbation:^ , Riegber uses twice the period o£ 
the Ci'ualyzed) arc as fundamental, and all necessary subharmonics 
of this to describe tlie variations of 'Kepler* element combin- 
ations. Using tiiese (arbitrary) frequencies, he calculates 
(by quadrature) theoretical amplitudes for the element combina- 
tions in the style of a Courier analysis. The advantage of 
the method is that it (apparently) separates effects perfectly; 
only orthogonal functions are determined. The disadvantage 
is that tlie natural frequencies are not all simply related to 
I even) a well chosen arc length. Therefore, a large number of 
terms may have to be evaluated to resolve a variation of a 
few close frequencies. More serious may be the i-estriction that 
only the amplitudes are determined by the method. Geopotential 
phase information is lost in defining two boundary values for 
each element from the orbit data. Nevertlieless , impressiv’e 
results have been achieved with this method (Riegber, 1974,1975) 
and an agreement with our analysis can be demonstrated. 

Riegber (1973) analysed optical data in a GBOS-II - 

arc (S.S days long) for the amplitude of the resonant 

variation oi w+M+^‘i!. The measured amplitude is related to a 

condition equation (a calculated quantity) involving all 
^ t ll 

13 order harmoni.cs. Unfortunately, no direct check can be 
made with our influence coefficients since ours are for 
sine and cosine terms independently. They contain phase as 
well as amplitude information. However, the amplitudes of 

1 T 1 / ^ 

our lumped coefficients ICC*) "+ (S*) "J can be compared to 
Riegber ' s . 


In Table 14 we make this comparison for 3 fields: 

1. Our GEOS -11 field from individual arcs (Table 7, 
along track) 

2. The SAC SE2 
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TABLb 14 COMPARISON or ':i.uS-2 RESULTS WITH 
A RIEGBER CONSTRAINT 


FOURIER' 

PERIOD 

(DAYS! u)® 


DATA ARC IS 6.8 DAYS LONG 


OBSERVED 

AMPLITUDE 

(10~^ RADIANS) 


CALCULATED AMPLITUDE (10~^ RADIANS) 

RIEGBER 

OUR INDIVIDUAL ADJUSTED 
SAP SE 2 ARC SOLUTION SOLUTION 


6.8 


261.2 300 4 16 303 303 


327 


CALCULATED AMPLITUDES ICVSM 
ANOIRATiO) ICVSM /CALCULATED 
RIEGBER AMPLITUDE 


66.95 66 91 62 32 

(6.4) (5.4) (6.3) 
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3 . 


Riegber's (^1973) corrected with IL£1^ (17,13)^, j. ■ 
(60 , -4 . 5) 

It is noted that our field agrees with Riegber's observation 
as well as liis whicli was 'fit* to tliis data. 

The computed lumped coefficients for u+M+n by our 
analysis are very close to those for the along track since it 
is dominated by t»i+M which is the same for botJi, But the 
correspondence of the amplitudes of our lumped values (using 
the various fields) with those computed from Riegber's 
condition eciuations is close but not exact. It is gratifying 
that the comparison is good. 
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' j 

Verification o£ Results Usint* GHOS-Il Dat a ; 

' I 

Another obvious means for assessing the accuracy of the ; 

derived values for the resonance coefficients found in '! 

i ( 

this analysis is to use these coefficients and compare orbit j 

determinations with them and other sets of coefficients. ( 

Five epochs were selected having a wide range of argument of ]\ 

perigee. The data reductions were repeated three different 
ways : 

• a solution was performed using the GEMO (Lerch, j 

1974) Gravity Model without any m=13 terms. This || 

was used as a basis for assessing the total 1| 

resonance information in the arcs. 

I 

• a solution was performed using the GEM6 Grav'ity || 

Model complete to (22,14). Three pairs of reso- || 

nance coefficients were adjusted - (14,13),, I 

L. , o 

(15,13)p Q and (10,13)., c. These solutions were ' 

used with those above to gauge the total resonance I 

j( 

information in the arcs. i! 

I 

y. 

• and lastly, the same arc was reduced using GEMb I 

ii 

without its 111=^13 coefficients, but added to 
GEM6 were the resonance coefficients recovered from J 

the multi-arc analysis (presented in Table 10). |l 

When this solution is compared to the one immediately jj 

above, one can get a fairly accurate measure of | 

resonance modeling obtained using the constraints || 

derived in this analysis. f 

[i 

I 

Table 15 summarises these results. I 


i! 

i 
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It is Idealized that this form of verification has cer- 
tain limitations. When one introduces six additional degrees 
of freedom to any data reduction, other errors *itag, 

solar radiation pressure, low degree and order geopotentialj 
will be partial!/ accommodated. Therefore, the solution using 
the complete GEM6 with three pairs of coefficients adjusting 
may yield results which include accommodation to these other 
error sources. Nevertheless, the level of resonance modeling 
obtained from our analysis can be inferred (pessimistically) 
from the results presented in Table 15, Our global solution 
models all but about 1.7*i of the IStli order resonance informa- 
tion in the GEOS -11 orbit. 

Table 15 uses weighted RMS as its measure of tlte quality 
of fit to the data. The actual data weights employed were: 

camera observations : 2*' 

range riaser) observations: 5m 

The laser range data was also sampled so that only 5U pts/pass 
was used iii tliis analysis. 
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TABLE 15. 

RESONANCE VALIDATION USING GEOS-II DATA 


RMS OF FIT TO THE DATA 


EPOCH 

OF 

ARC 

ARG 

OF 

PERIGEE 

(EPOCH) 

GEM 6 
W/OUT 
ANY 0 

lm«13) 

RES TERMS 

GEM 6 
ADJ. 

(14. 13) 0 
(15. 13) 

(16. 13) 

GEM 6 
USING 

DERIVED 0 
(m=13) 

RES TERMS 

© 

TOTAL 

RESONANCE 

CONTRIBUTION 

© - © 

ESTIMATED 

DEGRADATION 

USING 

DERIVED 

VALUES 

© - ® 

PERCENT 

performance 

OF 

DERIVED 

MODEL 

© + ® 

680509 

337® 

10.616 

0.840 

0.876 

9.776 

036 

ao% 


680407 

30® 

9.224 

0.776 

0.823 

8 448 

047 

0.6% 


680907 

140® 

3.348 

0857 

0.885 

2.491 

028 

1.1% 


690207 

252® 

14.739 

0930 

1 360 

13.809 

430 

3.1% 


680611 

282® 

9666 

0.871 

1.026 

8.796 

.155 

1.8% 



4 ^ 

tn 


RMS 

PERFORMANCE 


1.7H 




CONCLUSIONS 


Gravitational constraint equations have been derived 
i-ro)ii GEOS^TI data and a detailed analysis of the shallow 
resonance problem. These equations follow from elementary 
perturbation theory and the along- track constraint derived 
from them accounts for all but about of the 13th order 
resonant information in the tracking data. The equations 
are also in good agreement with recent comprehensive gravity 
models of SAG, GSTC and from Doppler data only, which use 
substantial amounts of GEOS -I I data. 

The goal of the analysis has been met; to derive 
simple constraints which account for nearly all of the shallow 
resonant information in satellite orbits. Tlie method makes 
it feasible to reduce compreliensive field models to the lumped 
coefficients for orbits or orbital arcs used in their solutions. 
The proper combination and extension of this constraint data 
should be straightforward, but is a task left for the future. 
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